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1 Introduction 



The study of steady states of non-equilibrium systems has motivated a 
lot of works over the last decades [1-9]. It is now well established that 
non-equilibrium systems exhibit in general long-range correlations in their 
steady state [4,10-12]. 

One of the most studied examples of non-equilibrium system is the one 
dimensional exclusion process with open boundaries [4,10,13]. The system 
is a one dimensional lattice gas on a lattice of L sites. At any given time, 
each site (1 < i < L) is either empty or occupied by at most one particle 
and the system evolves according to the following rule: in the interior of 
the system (2 < i < L — 1), during each infinitesimal time interval dt, a 
particle attempts to jump to its right neighboring site with probability dt 
and to its left neighboring site with probability q dt. The jump is completed 
if the target site is empty, otherwise nothing happens. The parameter q 
represents a bias (i.e. the effect of an external field in the bulk). The 
boundary sites i — 1 and i = L are connected to reservoirs of particles 
and their dynamics is modified as follows: if site 1 is empty, it becomes 
occupied with probability a dt by a, particle from the left reservoir, and if 
it is occupied, the particle is removed with a probability 7 dt or attempts 
to jump to site 2 (succeeding if this site is empty) with probability dt. 
Similarly, if site L is occupied, the particle may either jump out of the 
system (into the right reservoir) with probability (3 dt or to site L — 1 with 
probability q dt, and if it is empty, it becomes occupied with probability 
5 dt. 

The rates a, (3, 7 and 8 at which particles are injected at sites 1 and 
L can be thought as the contact of the chain with a reservoir of particles 
at density p a at site 1 and with a reservoir at density p^ at site L. The 
reservoir densities p a and pb are related to a, (3, 7 and S by (see appendix) 
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Pa = (1-la) 

l-g-l3-S+ v /(/3-<5-l+ g )2 + 4^ 

Pb = w^) (Llb) 

For q = 1, the bulk dynamics is symmetric and the model is called the 
Symmetric Simple Exclusion Process (SSEP) [2,3]. In the steady state, 
there is a current of particles flowing from one reservoir to the other (when 
Pa Pb) and the steady state density profile is linear. 

For < q < 1, the bulk dynamics is asymmetric and the model is 
called the Asymmetric Simple Exclusion Process (ASEP) [14-16]. When 
the densities p a and pb vary, the system exhibits phase transitions, with 
different phases: a low density phase, a high density phase and a maximal 
current phase [16-19]. On a macroscopic scale, the steady state profile 
is constant except along the first order transition line p a = pb < 1/2 
separating the low and the high density phases. 

In the large L limit, the probability Pl({p(x)}) of observing a given 
macroscopic density profile p(x), < x < 1 can be expressed through the 
large deviation functional J-({p(x)}; p a , pi,) by 

Pl{{p(x)}) ~ e'Lr({ P (x)};p a , Pb ) ( L2 ) 

This large deviation functional T is an extension of the notion of free en- 
ergy to non-equilibrium systems [20-22] . 

One can think of a number of distinct definitions of Pl({p(x)}) which 
all lead to the same T in the large L limit. Here, by dividing the sys- 
tem into k boxes of size L\, L2,. ■ ■ , (with X)i=i Li = L), we define 
QLi,L 2 ,...,L k (Ni, N2, . . . ,Nk) as the probability of observing in the steady 
state Ni particles in the first box, N2 in the second, . . . , Nj. in the last 
box. Then if we identify P L ({p(x)}) with qL u Lz,...,L k {Ni, N 2 , . . . ,N k ), one 
expects that for large L, (jl.2|l holds when 1 C L,; < L and Ni is the integer 
part of Lip(xi) where Xi = Yl)=i it 
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In [20,21,23-26], the following exact expressions of J-({p(x)};p a ,p b ) 
were obtained: 

In the symmetric case, i.e. for q = 1, it was shown in [20,21] that 

J r SSEp({p{x),p a ,Pb}) = SUp/ dx<p(x)\0g — T + (1 - log ' 



F(x) JO 

+ log( 

,Pb - Pa 

where the sup is over all monotone functions F(x) which satisfy 

F(0)=p a , F(l)=p b . (1.4) 

The auxiliary function F(x) which achieves the sup is the monotone solu- 
tion of the nonlinear differential equation 



F(x) J v rv " 6 \\-F(x) 



with the boundary conditions (|1.4|) . 

In the asymmetric case (i.e. for q < 1), the expression of J^({p(x)}; p a , pb) 
is given [25, 26] 

• in the case p a > p b by 

•^ASEp({/o(a;)}; Pa, Pb) = -K(p a ,p b ) (1.6) 
+ sup/ ^p(^)log[p(x)(l-F( a; ))] + (l-p(a ; ))log[(l-p(a ; ))F(a ; )], 

F(x) JO 

where the sup is over all monotone non-increasing functions F(x) 
such that F(0) = p a and F(l) — p b and 

K{p aiPb )= sup log[p(l-p)], (1.7) 

Pb<P<Pa 

As shown in [26], the function F(x) which gives the sup in (Ol) is the 
derivative of the concave envelope of Jq[1 — p(x')]dx' , whenever this 
derivative belongs to ]pb, p a [, and it takes the value p a or p b otherwise. 
As a result, when F(x) differs from p a and p b , it is made up of a 
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succession of domains where F(x) = 1 — p(x) and of domains where 
F(x) is constant (as F(x) is decreasing, it cannot in general coincide 
with 1 — p{x) everywhere). In the domains where F{x) is constant 
(and differs from p a and pb) it satisfies a Maxwell construction rule: 
when F(x) = C for 0<t<x<u< 1, its value is determined by 

(u-t)C = J [l-p(x)]dx (1.8) 
• and in the case p a < Pb by 



^ASEp({p(a;)}; p a , p b )) = -K(p a , p b ) + (1.9) 

ni n t, \ / dx p( x ) 1o s M^X 1 _ P")] + p( x )) lo s K 1 ^ pO))/ 5 ®] 

+ jf p(z) log [p(a;)(l - p 6 )] + (1 - p(z)) log [(1 - p{x))p b ] | . 
where 

K(p a ,Pb) = min [logp a (l - p a ),logp b (l - p b )] , (1.10) 

The goal of the present paper is to reconcile the expression valid in 
the symmetric case (|1.3f) with those valid in the asymmetric case <|1.6I1.9[) 
by calculating the large deviation functional in a weak asymmetry regime, 
which interpolates between the two, where q — > 1 as L — > oo with q = 1 — -j- 
The SSEP and the ASEP appear therefore as limiting cases of the results 
obtained in the present paper. 

The outline of the paper is as follows: in section 2, we summarize our 
results by writing several equivalent expressions of J- in the weak asymme- 
try regime. In section 3 we give the details of our derivation. In section 4, 
we show how the SSEP and the ASEP expressions can be recovered in the 
limit A — > and A — > oo. The large A limit is somewhat reminiscent of the 
WKB method and the calculation of the position 1)1. 8JI of the plateaux in 
the Maxwell construction of the function F(x) have an origin very similar 
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to the Bohr-Sommerfeld rule in the WKB method [27,28]. In section 5 we 
extend the range of validity of our results and show in particular that they 
remain true when detailed balance is verified. 



2 Main results 



We consider here a weak asymmetry regime defined as a situation where 
q — > 1 as L — > oo, 

A 



5 = 1- 



(2.1) 



keeping A fixed. 

For technical reasons which will become clear at the end of section 3.3, 
our results are limited to the case 



A > 



and 



Pa > Pb 



(2.2) 



In section 5, we will discuss some extensions to a broader range of param- 
eters. 

Our main result is that in the weak asymmetry regime, the large devi- 
ation functional is given by 

HiP( x )h Pa, Pb) = -K x (pa,Pb)+irf \ y(0) log Pa +y(l) log 1 

\v(x)\ 1 — Pn 



dx 



1 - e~ Xy 

log + plogp + (1 - p) log (1 - p) 



Pb 



(2.3) 



+ (1 - p + y') log (1 - p + y') + (p - ?/) log (p - y') 
where the inf is over all continuous positive functions y(x) satisfying 

p(x) - 1 < y'(x) < p{x) . 
We will show at the end of section 3.5 that the constant K\{p a ,Pb) is given 

by 



K\(p a ,Pb) = l0g(J) 



dp 

Mi - P ) 



log 1 



j 



(2.4) 
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where the parameter J is solution of 

dp 



„, J-Xpii-P)- 1 ■ (2 - 5) 

The parameter J is in fact related to the steady state current j by (see 
section 3.5) 

J= lim Lj (2.6) 

L — >oo 

Expression (|2.3|) for T can be rewritten in a form which interpolates 
between the symmetric (|1.3f) and the asymmetric <|1.6|1 cases (see section 
E3> : 



J dx j/>log J; + (l-p)log:j!- 



F 



+ log (F(l - FU - F') + rlogf ; — y ) 1 (2.7) 

sv v ' ' XF(l-F) & \ F(l-F)X-F'Jj 1 ' 

where the function F(x) is the solution of the differential equation 

[F - p)F' 2 + F(l - F)F" + XF(1 - F)(F - 1 + p)F' = (2.8) 

with the boundary conditions 

F(0) = Pa F(l) = Pb (2.9) 

In the range of validity of our derivation (|2.2() this differential equation has 
a unique solution, and this solution is monotone (see section I8.4J1 . 
Actually, if we consider the right hand side of H2.7(l as a functional of 
function F, then (|2.8|l appears to be the condition that F maximizes this 
functional under the constraint l|2.9|l (see the end of section I3.4|) , leading 
to 



F{{p{x)},Pa,Pb) = SUp 



-K\(p a ,p b ) 



dx^p\Og^ + (1 -p)l0gy 



F 



l-p_ 

-F 



log (F(l - F)X - F') + — — log 



XF(1 — F) b \ F(1-F)X-F' 



(2.10) 
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where the sup is over all decreasing functions F satisfying 1|2.9|) . 

As J 7 is a sup over convex functions of p, it is a convex function of p in 
domain (12.21) . 

A by-product of 12. 7|) is (see section I3.5[l that the most likely profile 
p(x) is solution of 

Depending on the boundary conditions this leads either to a tan profile, a 
tanh profile, or a coth profile. 



3 Derivation 

3.1 The matrix method 

The equal time steady state properties of the ASEP can be exactly calcu- 
lated using the so-called matrix method [14]. Let us consider a microscopic 
configuration defined by its occupation number {t{\ where t% = 1 when 
site i is occupied by a particle, and otherwise. It can be shown that the 
steady state probability of such a configuration for a lattice of L sites can 
be written as 

with Zi,(q) being a normalization factor defined by 

Z L (q) = (W\(T> + -E) L \V) (3.2) 

where D and E are two operators fulfilling the followings algebraic rules: 

DE - gED = D + E , (3.3a) 
{PD-SE}\V) = \V), (3.3b) 
(W|{aE- 7 D} = (W\ . (3.3c) 
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These rules <|3.3a|) - (|3.3c|) allow the computation of all equal time steady 
state properties without the need of finding an explicit representation. 

The two point correlation function < TiTj > (where the symbol < . > 
stands for the average with respect to the steady state probability) is given 

by 

(W\(D + E) i ~ 1 D(D + Ey- l - l D{D + E) L ~i\V) 



Z L (q) 

and the steady state current j between site i and i + 1 is given by: 



(3.4) 



j =< Ti(l - n+i) > —q < (1 - n)n + i > (3.5) 

Using expression l|3.4|l for the correlation function and the algebra rule 
H3.3afl . one gets 

] - Z L (q) • (3 - 6) 
Clearly the current does not depend on the site i, as it should, due to the 
conservation of the number of particles. 

If we divide the system of size L in k boxes of size L\, L2, ■ ■ ■ Lk the 
probability of finding N\ particles in the first box, N2 in the second, . . . and 
Nk in the last box is given by 

{W\X Ll {N x )X L2 (N 2 ) ■ ■■X Lk (N k )\V) 



qL u L a ,...L h (Ni,N 2 ,~- ,N k ) 



Z L (q) 

(3.7) 



where Xi (N) is the sum over all products of / matrices containing exactly 
N matrices D and / — N matrices E. 



3.2 A representation for D and E 

All physical quantities such as (|3.5|l or l|3.7|l do not depend on the 

representation of the matrices D and E and of the vectors \V) and (W\ 
which satisfy the algebra H3.3afl - H3.3cfl . Several representations have been 



9 



used to solve (|3.3al3.3c[) [14-17,29]. We choose in this section a particular 
representation which will be convenient for the remaining of our derivation. 
If we write the operators D and E as infinite matrices of the form 



D 



1-9 



E 



1-d 1-q 

1-dq l-q 2 

1 - dq 2 1 - q 3 

1-e ••• ' 

1-ed l-eq 

1 - edq 1 - eq 2 

'■. *•. '■. '•. 



(3.8a) 



(3.8b) 



we find that they satisfy the algebraic rule l|3.3a[l for arbitrary choices of d 
and e. 

Let call {|^)}n>i the vector of the associated basis. If we look for 
vectors \V) and (W\ of the form 



1 - Pa 
Pa 



71=1 

n— 1 ^ 



PbJ {q;q)n- 



(3.9a) 
(3.9b) 



where p a and p^ are for the moment arbitrary and the symbol (x; q)i stands 
for the g-shifted factorial defined by (x; q)o — 1 and 



xq 



for k > 



(3.10a) 



fc=0 



one can check that (|3.9a|) and l|3.9b|l fulfill the algebraic rules (|3.3b|l and 
l|3.3c|) if the parameters p a , Pb, d and e satisfy: 



Pa 1 - Pa 



= 1-9 



1 — Pb Pb 



-~ = l-q 



(3.11a) 
(3.11b) 
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d 



S(l - Pb) 
Ppb 



(3.11c) 



and 



IPa 



(3. lid) 



a{\ - Pa) 



Note that the parameters p a and pt, defined by Ijl.lafl and Ijl.lbfl . that 
we interpreted as the reservoir densities, are the unique solutions of the 
equations i|3.11a[l and i|3.11bjl such that < p a < 1 and < p& < 1 (this is 
why we use the same symbol in (3.9) and (1.1)). 

If we use this representation when q < 1 (which is not a restriction due 
to the left-right symmetry), we see that the condition for (W|X|V) to be 
finite when X is an arbitrary product of D and E is that p a > Pb, leading 
to (|2.2JI . In sectional we will show that this representation remains valid 
for some part of the domain q > 1, thus allowing us to extend the result of 
section 2 . 

3.3 The sum over paths and the derivation of (J2.3|) 

The basic idea of our derivation of (|2.3() is to expand the matrix products 
such as H3.1(l or Ij3.7|l as a sum over paths [12], in much the same way as 
the path integral formulation of quantum mechanics. 

Consider the set of discrete walks w of k steps. Let rii(w) (0 < i < k) be 
the integer position of the walk after the ith step. The walks we consider 
remain positive rii > and their increment at each step satisfies 



n t > 



and 



-Li < m - rij_i < L L 



(3.12) 



From l|3.7|l . we deduce 



1 



qL u L*,...L h {Ni,N 2 ,--- ,N k ) 



Z L (q) 



^(W\n )(n k \V)Y[(n i - 1 \X Li (N i )\n i ) 



W i — 1 



(3.13) 



where Xi(N) has been defined in (|3.7|) . 
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In the large L limit (with q — 1 — l|l.la|l and l|l.lbj) become 



Pa = — — 

a + 7 

and (|3.11d|) . jXTTcl) 



0(- 



d = i + 0(~) 



l 



e=l + 0(-) 



(3.14) 



(3.15) 



To compute (|3.13H when 1 -C Xj <C £, let us evaluate (n\X[(N)\n') when 
1 < I « i. When n, n' are of order L with |ra — n'\ < I <C L, we see that 
for m and m' in (min(n, n') — /, max(n, n') + 1) all the non-zero elements of 
matrix (m|D|m') and (m'|E|m) are equivalent to 



|D|m') - (m'|E|m) ~ 1- e 



A 



(3.16) 



The computation of (n\Xi(N)\n'} is thus reduced to an enumeration prob- 
lem. This leads to 

/ , , , W j_ N 

7J_ 



<n|Jr,(A0K> 



N 



£ 

nj r —n-—n' — n 

(3.17) 

where (^) is the number of words of length I with N matrices D and I — N 
matrices E, n+ is the number of matrix elements of the form {m\T)\m + 1) 
and n_ the number of matrix elements of the form (m|E|m — 1). 
Looking at the values of n + and n_ which dominate (|3.17() one obtains 

A 



\og(n\Xi(N)\n') = I 



log - + log (1 



-plogp- (1 -p) log(l -p) 



(1-P + y') log (1 - p + y') - (p - y') log (p - y') + o(l) (3.18) 



where y, y' and p are defined by 



V= L 



and 



N 

P = T 



(3.19) 



As -y 1 — » while Lj — > oo, one can associate to each walk w a continuous 



function y(x) 



»i(w) 



(3.20) 
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with Xi = X)j=i if- For each walk y(x) and each density profile p[x) such 
that 



y{x) > 



and 



p(x) - 1 < y'(a;) < p(z) 



(3.21) 



(these restrictions on the path y come from n.i(w) > (see (13.12(1 ) and the 
condition for (n\Xi(N)\n') ^ (see CSH)), let us define G({y(x)}, {p{x)}) 

by 



G({y(x)}, {p(x)}) = „(0) log + y(l) 



log 



1 - Pb 
Pb 



dx 



1 - e~ Xy 

log + plogp + (1 - p) log(l -p) 



+ (1 - p + y') log (1 - p + ?/) + (p - y') log (p - y') 



(3.22) 



Then expression (|3.13|l leads to (after replacing the sum over w by a sup 
over w) 



logq LuL2 ,.. Lk (Ni,N 2 r-' ,N k ) ~\ogP L ({p{x)}) 



isup [K x {p a ,p b ) - Q({y(x)},{p(x)})] 

V 



(3.23) 

*)})] 
(3.24) 



with K x (pa,Pb) given by 

K\(p a ,Pb) = lim logL 



logZ L (l-£) N 



(3.25) 



which is the expression l|2.3|l . 

In order to compute K\{p a , pb), we could use a similar method to esti- 
mate Zl (q) in (|3.25|) ; we will rather use the property that the most likely 
profile p must verify J 7 ({p(x)}; p a , Pb) — 0, so we delay the computation of 
K\{p a ,Pb) until the end of section 3.5. Of course both methods give the 
same result. 
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3.4 Derivation of (I^TH) 

As the functions — log(l — e - ^^), (p(x) — y'(x)) log (p(x) —y'{x)) and 
(1 — p{x) + y'(x)) log (1 — p(^) + y'(x)) are convex functions of y for every 
value of x , the function G({y(x)}, {p(x)}^j in 13.2211 is a sum of convex 
functions, and is thus convex. 

Therefore there is a unique walk y m (x) which minimizes Q ({y(x)}, {p(x)}^j , 
so 

T ({p(x)} ; Pa , Pb ) = g({y w (x)}, {p}) - K x ( Pa , Pb ) (3.26) 

This minimum is not reached at the boundary of (|3.21|) . Thus y m (x) is the 

. . , .. t 9G({y},{p}) 
unique stationary point, solution ot 



dy(x) 



= and it satisfies: 

9 = Sm 



^~y'^) -A g^L = (3.27) 



(1 - p(x) + y' m (x))(p(x) - y' m {x)) 1 - e - A »-(») 
with the boundary conditions 

1 - p(0) + j4(0) - p a 1 - p(l) + = Pb (3-28) 

in addition to the general condition (|3.21(1 . 

To obtain an expression similar to l|1.3|) . we rewrite l|3.27|l using the 
function F(x) defined by 

F(x) = l-p(x)+y' m (x) . (3.29) 

This lead to 

F'( T \ P -Ay m (x) 

= A , x (3-30) 



F{x)(l-F(x)) i _ e-A^Cx) 

with the boundary conditions 

F(0)=p a F(l)=/»6 (3-31) 

The conditions on the walk H3.21(l imply 

< F(x) < 1 (3.32) 
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and as y > (see (EHT|) and KFM ) 

F'(x) < (3.33) 
By eliminating y between <|3.29|) and (|3.3U|) we get (|2.8[) . 

The expression of the large deviation functional can be rewritten in 
terms of F instead of y m in (|2.3() . From l|3.30[l , we see that 

AF(1 - F) 



log(l-e- A ^) - log- 



(3.34) 



\F(1 — F) — F' 

If we integrate by part the term dxy'(x) log in (|2.3f) . we get 

l-p + y 111 



dxy log — 

P- V 



2/log- 



dx 



y(p' y") 



{I- p + y>)[p-y>) 
y(l)log-^-t/(0)lor 



(3.35) 



1 -P6 
1 F / 



log 



1 - Pa 

—F' 



7 (3-36) 



\F{1 -F) ° \F(l -F)-F 
and thus lf2~3j) leads to lETFl) . 

Let us now justify (|2.1()|1 . If we define 

H {{p}, {<P}) = -K x (p aiPb )+ f \p(x) \ogp(x) + (1 - p(x)) log (1 - p(x)) 



<t>'{x) 
A 



log 



-(l-p(x))(f)(x)-log(l + e^< 
(-*'(*)) + (l - ^) log (A -^'(s))} , (3 



37) 



expression (|2.7fl implies that when F is solution of (|2.8|) with condition 
(J2jU then 



•^(M; Pa, pb) = W (^{p}, jlog ^~ 



(3.38) 



As - log(l + e*< a >) and £M log ( - <t>'{x)) + (l - log (A - <f/(x)) 

are concave function of <fi for </>' < 0, H ({p}, {</>}) is a concave func- 
tion of (p. So 1|2.8[) is equivalent to the condition for F to maximize 
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H ({p}, | log (t^f) }) under the constraint (HJ and F' < 0, 

^({^{log^)}) _ 



OF 



So (|3.38ll can be written as 

F({p(x)};p a ,p b ) = sup 



H ({p}, {log 



(3.39) 



(3.40) 



where the sup is taken over decreasing functions F with the condition l|2.9|l , 
leading to (l2~lT)l) . 

3.5 The most likely profile 
From l|2.1U[l . we get 



dT{{ P {x)};p a ,p b ) _ p(l - F) 
dp(x) 8 (1 - p)F 

so that for the most likely profile p{x) 



(3.41) 



(3.42) 



F{x)=p{x) . 
Equation 1)2. 8|) thus becomes 

p" = Xp'(l-2p) (3.43) 

with boundary conditions 

p{0) = Pa and p(l)=Pb (3.44) 

Integrating (|3.43(l once, we get constant parameter J 

p = Xp(l -p)-J (3.45) 

The boundary conditions (|3.44|) determine J as given by equation l|2.5|l and 
pby J2HJ. 

One can show that (|2.7|l implies that 

1 



< TjTj + l > =< Ti X T i+ l > +0 



L 



(3.46) 
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(see for example [23] for the case A = 0). 

In particular < r^(l — Tj+i) >~< >< 1 — Tj+i > and (|3.5() for g = 1 — A 
leads to 



A 



i = tK 1 -p)-t + °( ^ 



(3.47) 



Comparing (|3.47|l to (|3.45(l leads to relation i|2.6[l between J and the current 
3- 

Depending on the reservoir densities p a and p^, one gets various expres- 
sion for the most likely profile p(x) 



• when p a - p b > A (p a - \) (pb - -j) then J > \: 



-, \ 1 /J- A/4 
p(.x) = J - y ^ tan 



y/\(J-\/±)(x - X ) 



(3.48) 



and when p a — pb < X (p a — §) (p& — ^) then J < ^ and 
I /A/4- J 



p(x) 



A 



■ coth 



VA(A/4- J) (a; - z ) 



(3.49) 



where Xq and J are chosen to satisfy 1(3.441) . 

Starting with the expression of the most likely profile, one can get an 
equation for K\{p a , p&) by writing that 



F({PY, Pa, Pb) = 



(3.50) 



Introducing equations 1(3.42(1 and ((3.45(1 in expression 1)2.711 and solving 
H3.50[> we get expression 1)2.4(1 for K\(p a , pb). 



4 Limiting cases 

In this section we show how previously known expressions (|1.3|l - l(1.7|l can 
be recovered as limiting cases of the results 1(2.312.11(1 of the present paper. 
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4.1 The SSEP limit 

Let us first consider the small A limit. Expression (|2.5[) for the current can 
be expanded in powers of A. 



. = Pa- Pb + A / Pa+p b _ pl+ papb+ P 2 b \ + Q ^ ( { ( . 

L L V 2 3 

in agreement (when A — > 0) with j = Pa ~ L Pb of [3]. 
The large deviation functional (|2.7|l becomes for small A 



T = 



£d*{plog(|)+(l-p)log(i-| 



F 

+ iog(-F') XF[ l pi F) | + 1 - K x ( Pa , Pb ) + 0(A 2 ) . (4.2) 

The leading order in A agrees with the SSEP expression (|1.3|) as the con- 
stant K\{p a ,pb) given by (|2.4|l becomes 

Kx{p a , Pb) = log(p a - pb) + 1 + 0(A) (4.3) 

and the equation 11.5fl for F is the limiting case of 1)2. 8|) when A — » 0. 
Furthermore one can check that the most likely profile (|2.11|) becomes linear 
as J = Lj — > p a — pb in the small A limit. Thus the results known for SSEP 
are recovered from our general A case in the limit A — > . 

4.2 The ASEP limit for p a > p b 

Let us now see how the strongly asymmetric case (|1.6f) can be recovered 
from (|2.7|l for large A. We see that for large A, the solution J of (|2.5|) is 
J ~ max A,o(l — p), so that 

Pb<P<Pa 

j - sup jp(l - p) when p a > p b (4.4) 

p b <p<p a 

in agreement with [18,30], and that K\(p a ,pb) ~ logJ in l|2.4(l so that 
reduces to lfO)| . 
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An interesting aspect of this large A limit is to see how the solution F 
of (|2.8|l becomes for large A the function F constructed from 1 — p{x) by 
the Maxwell construction explained after equation (|1.7|) . When A becomes 
large, (|2.8|) implies that F ~ 1 — p or F' ~ 0. Therefore one expects a 
succession of domains where F(x) ~ 1 — p(x) and domains where F(x) 
is constant, with the remaining contraints that F is monotone and p^ < 

F(x) < p a 

In a domain where F(x) ~ 1 — p(x), one can expand F in powers of -r 

l-2p p"" 



F = l-p+I 



O(^) ^ (4.r,i 



P (l - p) p' 

The condition F' < (|3.33|) implies that in such domains p' > 

In a domain t < x < u where F(x) ~ C is almost constant, neglecting 
F' 2 in O ( as < gives 



A / (C - 1 + p{x'))dx' 



(4.6) 



where B is constant over the whole domain t < x < u. 

The next question is to understand for large A the transitions between 
these different domains. 

A first possibility is that F(x) (which is monotone) becomes discontin- 
uous at such transition point. This is the case for example when 1 — p(x) 
is decreasing and discontinuous, implying that F(x) has a variation of the 
order of 1 (the discontinuity of p) over a range in x of order of j. This case 
can be analysed without difficulty, but we won't discuss it here. 

The other possibility is that for large A, F remains continuous but F' 
becomes discontinuous. This is what happens for example on figure ^ for 
large A there is a succession of domains where F — 1 — p{x) and where F 
is constant. 

Let us consider a domain where F(x) = l — p(t) = l — p[u) for t < x < u 
in the large A limit, surrounded by two domains where F(x) = 1 — p(x). 
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Figure 1: F(x) solution of i|2.8[l versus x for a given profile p(x) and reser- 
voir densities p a = 0.9, pb = 0.3 when A takes the value 1, 100 and 10000. 
The full-curve is 1 — p(x). The curves are obtained by minimizing (|2.3|) 
numerically for y{x) (discretized over 200 points). F(x) is then calculated 
from y{x) by relation (|3.29|) . For large A, we see a succession of domains 
where F = 1 — p(x) and where F is constant satisfying the Maxwell con- 
struction. 

At the transition point t, the solution of (|2.8|) takes a scaling form 

F(x) - 1 + p(x) = \[^G((x - t) yfJW) (4.7) 
where G is solution of 

G" = G(l - G') (4.8) 
For F to match the solution 1 — p(x) + 0(j) for x = t — and 1 — p(t) + o(l) 
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for x = t + 0, one needs that 

G(z) -> when z -> -oo (4.9) 

G(z) = z + o(l) when z -> oo (4.10) 

It can be shown that then G is solution of 

G' = l + wY-Ae - ^) (4.11) 

where W is the product logarithm function (also called the Lambert func- 
tion, see [31]) defined here as the largest real solution of 

W{x)e w{x) = x (4.12) 

Condition (|4.9II determines A = — e _1 (as W (— e _1 ) = —1). Using that 
for x small W(x) ~ x, the limit of G for large z can be computed 

G(z) -z~- er— dz' (4.13) 

It gives for the asymptotic regime (when « i - i « 1) of F' 

F> {x) „ _^) e -^(,- t ) 2 (4J4) 
e 

At the other boundary it of the domain, i* 1 has a similar scaling form. 

F(x)-l + p(x)=-^^G(-(x-u) s /p 7 {ujX) (4.15) 

with the same G solution of (|4.8I) . This gives for the asymptotic regime 
■j- « -(x - u) < 1 

* _^M e -^(-«) 2 (4.16) 
e 

For the asymptotics of (|4.(>|> to match with l|4.14J) as x — > t 

F'(x) - Be ^ (4.17) 

and with 14.16fl when x — > it 

F'(x) - Se[- A ^" (c-i+*'))^]-^(x-«) 2 ( 4 18 ) 
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one needs that to the leading order in A 

(C - 1 + p{x'))dx' = (4.19) 

which is the Maxwell construction (|1.8fl . We see that the constant C, 
the value of F(x) in a domain where F is constant, is determined by an 
expression Ij4.19|l which is obtained from two matching conditions at the 
boundaries of the domain. This is very similar to the Bohr Sommerfeld rule 
which determines the energy levels in the WKB method [27,28]. So the 
Maxwell construction here has a mathematical origin similar to the Bohr 
Sommerfeld rule. 



5 Extension of our results 
5.1 Extension of the domain (12.21) 

We are going to show that our results of section 2, initially derived in the 
domain (|2.2|1 remain valid for 

A>0 eX l- Pb Pa Kl (51) 

Pb 1 - Pa 

The representation of the algebra for D, E, \V) and (W\ introduced 
in section 3.2 remains valid for some range of the parameter p a and pb 
when q > 1. Indeed, for large n the non-zero matrix elements of the kind 
(n\T)\n') and (n'|E|n) behave then like q n . Furthermore, ~ (y)" 

for large n. Thus, using (|3~lTc] and (|3.11d|) . we get that for any product 
X of L matrices D or E (and any sum of such product) (W^ln) (n^ly) ~ 
l ?' L ~ 1 ) an d tne condition for to be finite is 

S^" 1 < 1 • (5-2) 
When q = 1 — ^ and in the large L limit, this leads to (using l|3.14fl ) 
A<0 e- x l ~ Pa Pb <1 (5.3) 

Pa 1 - Pb 
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So all the content of sections 3.3 to 3.6 remains valid, leading thus to 
formulas (|2.3|I - H2.11|I . The only change is that in l|2.10[l . the sup is now over 
decreasing functions F such that for every x, XF(x)(l — F(x)) — F'(x) > 0. 

In order to recover the expressions for A > 0, we use the left-right 
symmetry of the system, replacing A by —A, x by 1 — x, p a by pb, ... so 
that condition l|5.3|l becomes (|5.1|) . 

When condition l|5.1|l is fulfilled, J given by (|2.5|) is negative, and thus 
there is a current j going from the reservoir with the highest density pb to 
the reservoir with the lowest density p a , despite the external bias q. The 
most likely profile, solution of l|3.45|) and (|3.44(l is now given by 

p{x) = \ + 7 A/4 A ~ J tanh v/A(A/4- J)(x - x ) (5.4) 
5.2 The detailed balance case 

We show now that (|2.7|l remains also valid when the boundary parameters 
a, (3, 7 and S are such that detailed balance is satisfied. Detailed balance 
means that the probability of observing a transition from a microscopic 
configuration C to another C is equal to the probability of observing the 
reversed transition (from C to C). 

Let {ri} be the occupation numbers of a given microscopic configura- 
tion. Detailed balance corresponding to a jump of a particle between sites 
k and k + 1 means that 

P ({ti, • ■ • ,T fc _i,0, l,T fc+2 , . . . ,t l }) = q~ 1 P({T 1 , . . . ,r fc _i,l,0,r fe+2 , ■ ■ ■ ,r L }) 

(5.5) 

where P({tj}) is the steady state probability of configuration {Tj}. 
The detailed balance relation at the left boundary is 

P({1,T 2 ,...,T L })=-P({0,T 2 ,...,T L }) (5.6) 

and at the right boundary 

P({r 1) ...,r L _ 1) l}) = |p({ti,..., 7i_i,0}) (5.7) 
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Starting from a configuration with occupation number {Ti}, one can always 
use (|5.5|l and l|5.6l) to calculate the weights of all configurations by removing 
particles at the left boundary. 

p ({n}) = f[(-^y P({0,0,...,0» (5.8) 

i=X ^ ^ 

For general values of a, (3, 7 and 5, these weights do not satisfy 1)5.7(1 and 
are not steady state weights. If wc insist however that l|5.8|l satisfies also 
1(5. 7jl . we get 

which is the detailed balance condition, and if l|5.9|) is satisfied, we get 

P({n})= T y V \ (5.10) 



The average density at site i is thus given by 

1 

LI 1 + q 



P[t) = i , ( 5 - n ) 



In the weak asymmetry regime 1)2. ljl . expression 1)5.911 and 1)5.11)1 become 

1 (5.12) 



X 1 - Pb Pa 



Pb 1 - Pa 

and 

p(z) = x (5.13) 

whereas (15.1011 leads to the following expression for the large deviation 
functional 

F({p(x)}; Pa , Pb ) = J 1 jplog j|| + (1 - p(*)) log } <** • 

(5.14) 

Note that 1)5.12(1 corresponds to the boundary of the range of parameters 
1)5.1)1 where we have shown 1)2.7)1 to be valid. Although our derivation of 
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12.7[l was not a priori valid when detailed balance (|5.12|) is satisfied (see 
H5.ll) ). we are going to see now that (|5 . 131) and l|5.14|l can nevertheless be 
recovered. 

When detailed balance (|5.12|) holds, one can check that 

F (x) = P(*) = ~ h Pa S — x (5-15) 
p a + (1 - p a )e Ax 

is solution of lj2.8|l for arbitrary p(x). So Ij5.12|) implies that F(x) does not 
depend on p(x). As F(x) given by H5.15|) satisfies 

F' = \F(l-F) , (5.16) 

one can see that H2.7J) reduces to Ij5.14|l . 

On can also check that when detailed balance is satisfied the current (|2.5I2.6|) 
vanishes and that J = in 12.5( 1 is equivalent to l|5.12|) . 

6 Conclusion 

In the present work, we have obtained the expressions l|2. 31 12.7112. 1(J|) for the 
large deviation functional of the one dimensional simple exclusion process 
in the weak asymetry regime H2.ll) . Our analysis of the limiting cases 
(A — > and A — > oo) has shown that these new expressions are consistent 
with previously known expressions for the SSEP and the ASEP. 

For technical reasons, our derivation is limited to some ranges of pa- 
rameter (|2.2f) , 15.1(1 or 15.12(1 . It would of course be useful to know what 
happens in the other ranges of parameters, if our results remain valid or 
not and how the ASEP result ((1.911 can be recovered. 

The derivation of our results, based on the matrix representation of the 
steady state, uses strongly that the steady state weights of the configu- 
rations can be written as sums over paths of the weights of these paths. 
We used a similar idea recently to study the density fluctuations in the 
TASEP [12]. These paths have so far a purely mathematical origin and it 
would be of course interesting to give them a physical interpretation. 
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Another interesting question would be to see whether the results of 
the present paper could be understood using the macroscopic fluctuation 
theory [21,24]. 



A Appendix: Definition of the densities p 
(ll.lal) and ( jl.lbl) of the reservoir 



When the boundary parameters a, (3, 7 and 6 satisfy a certain relation 
f ljA.41) below), the steady state is a Bernoulli measure at density p and one 
can consider that the two reservoirs are at this same density p. 

When the steady state is a Bernoulli measure at density p a , the steady 
state current l|3.5f) in the bulk is given by 

3 = (I - q)p a (l - p a ) (A.l) 

and at the left boundary, by 

j = a(l - p a ) - 1Pa (A.2) 

The conservation of particles implies that (|A.1(1 and (|A.2(1 should coincide 
and this gives condition (|3.11a|l for p a . If one repeats the same argument 
at the right boundary, one gets that p b should satisfy 

(1 - q) Pb (l - Pb ) = I3p b - 6(1 - Pb ) (A.3) 

leading to equation l|3.11b(l . The solutions of (|3.11a[l and (|3.11b|) (satisfying 



< p a < 1, < p b < 1 ) are given in (11.1 aft and 1)1. lbl) . One recovers in 
particular p a = mm(j^,l), p b — max(l — j^,0) when 7 = S = as 
in [25], and p a = p b = when q — > 1 as in [23]. Comparing l|3.11a|) 
and p.llbp . one can check that for the two reservoirs to be at the same 
density p (i.e. for p a = p b = p), the boundary parameters should satisfy 

(a + (5)(/3 + 7 )(l -q) = (af3 - j6)(a + (3 + 7 + 5) . (A.4) 
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Let us now verify that when p a = Pb = P, the steady state measure is 
indeed a Bernouilli measure at density p. Consider a configuration {t{\ with 
the steady state probability P. The probability of leaving this configuration 
during the time interval dt is 

[(1 + 7)n + a(l - n) + n c (l + q) + 5(1 - t l ) + (J3 + q)T L ]Pdt (A.5) 

where n c is the number of clusters of particles in {t{\ which do not touch 
a boundary. 

If the steady state is Bernoulli at density p, the probability of entering the 
configuration {t^} is 



[n (q + a- P - J + 7 _^_(l-r 1 ) + n c (l + g ) 

V p / 1 -p 



+ r z (l + tfL-^) +;5 _£_(i_ Ti )]p^ (A.6) 

V p / 1 - p 

For (|A.5(1 and (|A.6|I to be equal for any {n}, one needs that 

1 + 7 — a = q + a - 7 — ^ — (A. 7) 

P 1 - P 

-* + /3 + g = 1 + 5^-/3-^- (A.8) 
P 1 - P 

a + 5 = 7— + /3— P — (A.9) 
1 - P 1 - P 

Comparing (|A.7|I and (|A.8|) to (|3.11a|) and (|3.11b|) . one sees that 

P = Pa= Pb (A.10) 

whereas (|A.9|I is equivalent to the difference of l|3.11a|) and l|3.11bj) so is 

automatically satisfied. 
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